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Quantum anomalies give rise to new transport phenomena. In particular a magnetic field can induce an 
anomalous current via the chiral magnetic effect and a vortex in the relativistic fluid can also induce a current via 
the chiral vortical effect. The related transport coefficients can be calculated via Kubo formulas. We evaluate the 
Kubo formula for the anomalous vortical conductivity at weak coupling and show that it receives contributions 
proportional to the gravitational anomaly coefficient. The gravitational anomaly gives rise to an anomalous 
vortical effect even for an uncharged fluid. 
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Introduction. The effects quantum anomalies have on the 
hydrodynamics of otherwise conserved currents have recently 
received much attention. Two such effects are known: an 
anomalous magnetic effect LI] and an anomalous vortical ef- 
fect 12. 01- In the first an (external) magnetic field induces a 
current via the axial anomaly whereas the second is the gener- 
ation of a current due to the presence of a vortex in the charged 
relativistic fluid. These effects have been argued to lead to 
observable event by event parity violation and a charge sep- 
aration effect in non-central heavy ion collision at RHIC and 
LHC |41]. In the hydrodynamic constitutive relations these ef- 
fects lead to the existence of a new class of transport coeffi- 
cients. 

A first principle calculation of transport coefficients is pos- 
sible via so called Kubo formulas. The Kubo formula for the 
anomalous magnetic conductivity has been derived and ap- 
plied in isll whereas the one for the anomalous vortical con- 
ductivity has been established recently in [6]. They are 

(^AB = }^<^T,£ijn:^{JAJB)\oy=o, (1) 

""A = \imye,,„-^{J\T'^')U=^, (2) 

where 7^ are the (anomalous) currents and T''^ is the energy 
momentum tensor These Kubo formulas allow the calculation 
of the transport coefficients in the constitutive relations 

J\ = aS^e-j'djA'.+a^e-^'djVk, (3) 

where Af are the spatial components of a collection of gauge 
fields and is the local fluid velocity. The fluid velocity is 
taken here in a laboratory frame, i.e. we take vj. independent 
of the external magnetic fields. There will also be a non-zero 
energy flow which could be compensated for by going to the 
Landau frame. Since we are interested only in the total current 
generated by the external magnetic fields we will not do so. 
See also |01 for a discussion. We note that we can substitute 
the fluid velocities by the gravitomagnetic potential. To do so 
we go to the rest frame of the fluid defined by uf^ — (1,0,0,0) 
but switch on a gravitomagnetic field in the metric according 
to 

ds^ = dt^ + lAgdxdt - dj? . (4) 



Using this metric we can compute the local fluid velocity — 
{\,Ag) such that v = Ag. All these expressions are valid only 
up to first order in the external fields Af,Vk. Plugging this 
into the constitutive relation ^ and noting that A'^ sources 
r"' leads to the Kubo formula for the vortical conductivity. 

We will now evaluate the Kubo formulas ([T]), (|2|i for a the- 
ory of free right-handed fermions transforming under a 
global symmetry group G generated by matrices {TaY g. We 
denote the generators in the Cartan subalgebra by H^- Chem- 
ical potentials ij.a can be switched on only in the Cartan sub- 
algebra. Furthermore the presence of the chemical potentials 
breaks the group G to a subgroup G. Only the currents that 
lie in the unbroken subgroup are conserved (up to anomalies) 
and participate in the hydrodynamics. The chemical potential 
for the fermion is given by jx^ = Y.A IaI^a, where we write 

the Cartan generator = 1^^^ g in terms of its eigenvalues, 

f ^ ^ 

the charges q^. The unbroken symmetry group G is generated 

by the matrices g fulfilling 

Tlgii^^llfjlg. (5) 

There is no summation over indices in the last expression. 
From now on we will assume that all currents Ja lie in di- 
rections indicated in (|5]l. We define the chemical potential 
through boundary conditions on the fermion fields around the 
thermal circle fk, 'P^(t) = -e'^^^>I'^(T - jS) with j3 = 1/r. 
Therefore the eigenvalues of 5^ are i(bn + /i^ for the fermion 
species / with &>„ = nT{2n + 1) the fermionic Matsubara fre- 
quencies. A convenient way of expressing the currents is in 
terms of Dirac fermions and writing 

4 - E T^f^gyy+^f, (6) 

/.s=i 

/=i 

where we used the chiral projector y± = j{l ± 75). The 
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fermion propagator is 



1 



- tEa 



(8) 

(9) 



with /d)-^ = ia)„ + 11^, = {l,tq), q — and Eg = \q\. 
We can easily include left-handed fermions as well. They 
contribute in all our calculations in the same way as the right 
handed ones up to a relative minus sign. 

Evaluation of Kubo formulas. We will address in detail 
the computation of the vortical conductivities Eq. (|2]i and 
sketch only the calculation of the magnetic conductivities 
since the latter one is a trivial extension of the calculation of 
the chiral magnetic conductivity in ^ . 

Vortical conductivity: The vortical conductivity is de- 
fined from the retarded correlation function of the current 
Ja (■"■) ® ' ^^'^ '■^^ energy momentum tensor or energy current 
T°j{x') 0, i.e. 

Gjix-x') = ieij„ie{t-t')mx),T'jix%. (10) 
Going to Fourier space, one can evaluate this quantity as 



Gjik) 



•5 ^ijn tl" 



sffiq)y 



xS-f f{q + k) {fq^ + yHm^) 



(11) 



which corresponds to the one loop diagram of Fig. [T] 

The vertex of the two quarks with the graviton is ^ 5^ and 
therefore we find only contributions from the diagonal part of 
the group G. Our metric is g^v = diag(l, —1, — 1, —1). We 
can split GJ into two contributions 



Gl{k)^Gl,,Jk) + Gl,.,m, 



(12) 



which correspond to the terms 'f'q^ and jHS)^ in Eq. ( fTTT i re- 
spectively. We will focus first on the computation of GJ 



The integrand of Eq. ( fTTT i for GJ^^q.^ can be written as 



(OjT 



jqJ eiJ„tr[Y^,Yrvfy^ 



(13) 



r.u=± 



xA,{ia)^,q)Auiia)-f +iCO„,q + k)q^{q + k)^^ . 



From a computation of the Dirac trace in Eq. (fTjt one has 
two contributions 

eijntrlr^YYyfja^b'' =4Sij„{a'b'' + a%-), (14) 
Sijntr[Y^YYvfY5]a^'b''=4i{ajb„~a„bj) . (15) 



q + k 




Figure 1: 1 loop diagram contributing to the vortical conductivity 
Eq. lO. 



Using Eqs. ( fTSl ), ( fT4] i and ( fTsT i one can express GJ^q^j(A;) as 



d q 



r,M=± 



t — 



k'+q' 



tu 



-i- 



■ [q jk„ qnkj) 



X Ai{ia)^,k)Au{i(i)^ + ia)„,q + k). 



(16) 



At this point one can make a few simplifications. Note that 
due to the antisymmetric tensor e,y„, the two terms propor- 
tional to q' inside the bracket in Eq. (fTSI l vanish. Regarding 
the term £jj„q^k', it leads to a contribution ^ £ij„yk' after in- 
tegration in d^q, which is zero. Then the only term which 
remains is the one not involving e,j„. We can now perform the 
sum over fermionic Matsubara frequencies. One has 



^Y^A,{i(i}f ,q)Au{iG)f + i(0„,q + p) = 

tn{Eg - tjif) - un{Ey+k - m^) + \{u-t) 
i(0„ + tEg - uEg+k 



(17) 



where n{x) — 1/ (e^'- + 1 ) is the Fermi-Dirac distribution func- 
tion. In Eq. dTTl) we have considered that co,, ~ InTn is a 
bosonic Matsubara frequency. This result is also obtained 
in Ref. After doing the analytic continuation, which 

amounts to replacing i(0„ by ^0 + '£ in Eq- ( flTl l. one gets 



^1(0;) W 

LI 

tM=± 



A f 



d^q q^k„~{q-k)q„ 



, J {iKy EqE,i+k 

i{Eg - tji-f) - tn{Eg+k - upif ) + \{t-u 
kfj + i£ + tEg ~ uEg^i^ 



(18) 



The term proportional to '-^ j{t — u) corresponds to the vac- 
uum contribution, and it is ultraviolet divergent. By removing 
this term the finite temperature and chemical potential behav- 
ior is not affected, and the result becomes ultraviolet finite 
because the Fermi-Dirac distribution function exponentially 
suppresses high momenta. By making both the change of vari- 
able q -> —q~k and the interchange —t and f ^ — m in the 
part of the integrand involving the term ~tn{Eg+ii — uii-f), one 
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can express the vacuum subtracted contribution of Eq. ([18) as 



{q-kf 



n{E^-^-f)+n{E^ + nf) 



tt^± ko + ie + tE^ 



- uE, 



q+k 



(19) 



The result has to be proportional to k„, so to reach this expres- 
sion we have replaced q„ by {q-k)k„/l? in Eq. ( fTSl l. At this 
point one can perform the sum over u and the integration over 
angles. Then one gets the final result 



I6n^ k^ 



f=i 



'k^+4q{q + tko)] log 



dqqf^{q) 

-{q+kf 
a}-{q-kY 





a} 



1 



(20) 



where £1, = + ' £ + tEq , and (q) = n{Eg - ji^) + n{Eg + 
jX^) . The steps to compute GJ ^^jq^ in Eq. (fTSl l are similar In 
this case the Dirac trace leads to a different tensor structure, in 
which the only contribution comes from the trace involving 75, 
i.e. 



£ijr, tr [r^LYYvY^Y5]a''b'' = ii{a„ba - aob„ 



(21) 



The sum over fermionic Matsubara frequencies involves an 
extra io)^. Following the same procedure as explained above, 
the vacuum subtracted contribution writes 



kn ^r 



N ^ 



dq'£f;^{q,ko) (22) 
,t± 



4tqkko - (r - kl) {2q + tko) log 



f + kf 



{q-kf 



where f^{q,ko) = qf^{q) +tkQn(Eq + tjxf) . The result for 
G^{k) writes as a sum of Eqs. (|20] | and (|22] |. according to 
Eq. (fT2] i. From these expressions one can compute the zero 
frequency, zero momentum, limit. Since 



lim 



Hm E log ( 

the relevant integrals are 

dqqf^iq)^ I dqf^{q,ko:==0)^ 
Jo 



2k 

q 



(23) 



— T^. 
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(24) 



Finally it follows from Eqs. ( |20] i and (l22T i that the zero fre- 
quency, zero momentum, vortical conductivity writes 



CT4 = 



1 ^ f 



f=i 



(25) 



16;r2 



2n^ 



J^tr {Ta{Hb,Hc})^b^+-^T^\^{Ta) 



Both GJjq^.j and GJ ^^-^ lead to the same contribution in . 
Eq. ( |25] | constitutes our main result. The term involving the 
chemical potentials is induced by the chiral anomaly. More 
interesting is the term ^ which is proportional to the grav- 
itational anomaly |8] as we will show in the discussion sec- 
tion. The Matsubara frequencies &)„ — nT{2n + 1) generate 
a dependence on nT in the final result as compared to the 
chemical potentials, and then no factors of K show up for the 
term in Eq. ( |25] l. Left handed fermions contribute in 

the same way but with a relative minus sign. Left and right 
handed fermions do not mix. 

If instead of having taken the zero momentum limit at zero 
frequency, one took the zero frequency limit at zero momen- 
tum, the result would be 1 /3 of the result quoted in Eq. ( |25l l. 
The same factor appears in the magnetic conductivity when 
one interchanges the two limits ijstl • 

Magnetic conductivity: The magnetic conductivity in the 
case of a vector and an axial t/(l) symmetry was computed 
at weak coupling in |5]. Following the same method, we 
have computed it for the unbroken (non-abelian) symmetry 
group G. The relevant Green function is 



- 2^7-^^««^|7y {2K 
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X tr 



Sffiq)ySff{q + k)YJ 



(26) 



The evaluation of this expression is exactly as in f3\ so we 
skip the details. The result is 



87r2 



Y^tv{TA{TB,Hc})^ic- 
(27) 



In the second equality of Eq. ( l27T i we have made use of 
Eq. (|5]l. No contribution proportional to the gravitational 
anomaly coefficient is found in this case. 

Discussion. In vacuum the anomaly appears on the level of 
three point functions. In the presence of external sources for 
the energy momentum tensor and the currents this is conve- 
niently expressed through [8J 



bA 



(28) 



The anomaly coefficients are defined by 



1 



dABC = -[tT{TA{TB,Tc})R-tv{TA{TB,Tc})L], (29) 



bA = tl-{TA)R-tl{TA)L, 



(30) 



where the subscripts R,L stand for the contributions of right- 
handed and left-handed fermions. 

We have computed the magnetic and vortical conductivity 
at weak coupling and we find contributions that are propor- 
tional to the anomaly coefficients ( |29] l and (|30] |. Non-zero 
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values of these coefficients are a necessary and sufficient con- 
dition for the presence of anomahes 1 10]. Therefore the non- 
vanishing values of the transport coefficients ([TJ and ([2) have 
to be attributed to the presence of chiral and gravitational 
anomalies. 

Since the gravitational anomaly is fourth order in deriva- 
tives it is a bit surprising to find it contributing to first order 
transport coefficients. One possible intuitive explanation one 
could think of is that the gravitational field in the presence of 
matter gives rise to a fluid velocity m'' e.g. through frame drag- 
ging effects and that this might effectively reduce the number 
of derivatives that enter in the hydrodynamic expansion. 

The holographic calculation in AdS/CFT |6] did not show 
a contribution proportional to T^. This not surprising since 
only a holographic gauge Chem-Simons term was included. 
Holographic modeling of the gravitational anomaly calls how- 
ever also for inclusion of a mixed gauge-gravitational Chern- 
Simons term of the form AARAR^. 

We find a non-vanishing vortical conductivity proportional 
to ~ even in an uncharged fluid. In 111 ill similar terms 
in the vortical conductivities have been argued for as undeter- 
mined integration constants without any relation to the gravi- 
tational anomaly. 

It is also interesting to specialize our results to the case 
of one vector and one axial current with chemical potentials 
liR ^ II + Ha, llL ^ IJ- - IJ-A, chai-ges = ( 1 , 1 ) and q^^ = 
(1,-1) for one right-handed and one left-handed fermion. We 
find (for a vector magnetic field) 



Ma 
' 27r2 

27r2 



■"AV 



JL 

' 27^2' 
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n 



(31) 



Here Gyy is the chiral magnetic conductivity 15|], cr^ de- 
scribes the generation of an axial current due to a vector mag- 
netic field |J2,1 . Gy is the vector vortical conductivity, is 
the axial vortical conductivity and the only one sensitive to the 
gravitational anomaly. 

In L 1 3il enhanced production of high spin hadrons (espe- 
cially baryons) perpendicular to the reaction plane in 
heavy ion collisions has been proposed as an observational 
signature for the chiral separation effect. Three sources of 
chiral separation have been identified: the anomaly in vac- 
uum, the magnetic and the vortical conductivities of the axial 
current J^. Of these the contribution of the vortical effect was 
judged to be subleading by a relative factor of lO^'*. The 
term in dSTl i leads however to a significant enhancement. If 
we take ji to be the baryon chemical potential « 10 MeV, 
neglect jiA as in 11311 and take a typical RHIC temperature of 
T = 350 MeV, we see that the temperature enhances the axial 
chiral vortical conductivity by a factor of the order of 10"*. We 



expect the enhancement at the LHC to be even higher due to 
the higher temperature. 

Beyond applications to heavy ion collisions leading to 
charge and chiral separation effects lll3[ 11411 it is tempting to 
speculate that the new terms in the chiral vortical conductivity 
might play a role in the early universe. Indeed it has been 
suggested before that the gravitational anomaly might give 
rise to Lepton number generation, e.g. in fis'l. The lepton 
number separation due to the gravitational anomaly could 
contribute to generate regions with non-vanishing lepton 
number. 
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